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The copper-oxide based high temperature superconductors have complex phase diagrams with mul- 
tiple ordered phases. It even appears that the highest superconducting transition temperatures for 
certain cuprates are found in samples which display simultaneous onset of magnetism and supercon- 
ductivity. We show here how the thermodynamics of fluid mixtures - a touchstone for chemistry as 
well as hard and soft condensed matter physics - accounts for this startling observation, as well as 
many other properties of the cuprates in the vicinity of the instability towards "striped" magnetism. 
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The phase diagrams of conventional superconductors 
are usually simple, with no ordered phases competing 
with the superconducting state. By contrast, the high 
temperature superconductors have a number of com- 
peting phases which appear as the temperature is low- 
ered. One of the most astonishing manifestations of 
this competition occurs in the LaCuO family of ma- 
terials where coexisting magnetism and high temper- 
ature superconductivity has been reported. In 
Lai.6-xNdo.4Sr x Cu04, long period magnetic ("stripe") 
order (as detected by neutron diffraction) sets in at a 
higher temperature than the superconducting T c (and in- 
deed, charge "stripe" order appears at a still higher tem- 
perature |j]||]). In La2Cu04.i2 and Lai.ssSro.isCuO^, 
superconductivity and long-period magnetism appear to 
have the same onset temperature in the same bulk crys- 
tal! However, while muon spin relaxation data tells a 
grossly similar story Jl|| in Lai.6_ a; Ndo.4Sr a; Cu04, the 
corresponding data |7| for La2Cu04.i2 show that the 
magnetism is peculiar in that it resides in only a frac- 
tion of the sample, and its temperature evolution is due 
to growth of the magnetic fraction rather than an in- 
crease in the order within the magnetic fraction. We 
show how these observations can be understood from the 
classical thermodynamics of two-phase mixtures, which 
is applicable because of the well-documented tendency of 
antiferromagnets to expel holes 10-13]. 



Phase Diagrams with Competing Orders: The 

interplay between "stripe" magnetism and superconduc- 
tivity can be understood most simply by treating the 
liquids of mobile charge carriers in the high-temperature 
superconductors as fluids with a variety of ground states. 
As for other complex fluids, the coupling between the or- 
der parameters can lead to phases with mesoscopic den- 
sity modulations as well as diverse combinations of the 
order parameters themselves. To make this phenomenon 
explicit, we follow a standard paradigm of statistical 
physics and consider jl3| the simplest Landau free en- 
ergy, F, for two coupled order parameters, S and A, 
which represent the long period antiferromagnet and the 



superconducting order respectively. Spin rotation invari- 
ance and gauge invariance (A is a complex number whose 
phase cannot influence the free energy) imply that F is 
a function of \S ■ S\ and | A| 2 : 

F = F (ji,T) + a(p,T)\S-S\+(3(^T)\S-S\ 2 (1) 
+ a ( M ,T)|A| 2 + 6( M ,r)|A| 4 + 7 ( Ai ,T)|A| 2 |5.5|... 

where T is the temperature, \i is the chemical poten- 
tial for doped holes, ... represents higher order terms in 
powers of the ordering fields, and the various coefficients 
embody the effects of all the short-distance physics. The 
phase diagram is then determined (at mean-field level) 
by minimizing F with respect to S and A. While such 
a mean-field description ignores important fluctuations, 
especially given the fact that the high temperature su- 
perconductors are quasi- two dimensional, it provides a 
valid zeroth order way to examine the global structure of 
the phase diagram. 

From the macroscopic viewpoint adopted in the 
present paper, the parameters which enter the Landau 
free energy in Eq. (1) are purely phenomenological. How- 
ever, some insight concerning the microscopic physics 
can be inferred from the behavior of these parameters. 
In particular, if 7 > 0, superconductivity and long pe- 
riod magnetism compete, while if 7 < 0, they enhance 
each other. Indeed, Tranquada and collaborators || have 
concluded that static magnetism and superconductivity 
compete, and this is certainly intuitively sensible. Recent 
experiments on the behavior of vortex cores []l4f and on 
superconductivity-induced changes in the magnetic sus- 
ceptibility p5| of optimally doped La2-xSr a; Cu04 con- 
firm that 7 > 0, in agreement with these arguments. 
One main purpose of this paper is to show how even if 
7 > 0, magnetism and superconductivity can set in at 
the same temperature in a single sample. 

An important subtlety arises from the fact that, as in 
many experiments on classical fluids, it is the total num- 
ber of constituents of the fluid rather than the chemical 
potential that is fixed. For the cuprates, the constituents 
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are the charge carriers (doped holes) and their number 
is fixed by the chemical composition of the compound 
under study (e.g. the x in I^-^Sr^CuC^). Therefore /i 
must be determined from the implicit relation 

-x= dF/d^ = F^ + a'\S ■ S\+ a'\A\ 2 + ... (2) 

where x is the concentration of doped holes, / denotes 
differentiation with respect to /i, and S and A are the 
equilibrium values of the ordering fields as a function of 
T and fi. However, where two-phase coexistence occurs in 
the phase diagram there are values of fi at which dF/d/i 
has a discontinuity. In this case, Eq. (||) has solutions 
for fixed [i at two different values of x, x\ and X2- The 
equilibrium state for fixed x in this range consists of a 
two phase mixture, with the volume fractions of the hole 
rich (x — X2) and hole poor (x = x\) phases determined 
by the classical lever rule, f± — (X2 — x)/(x2 ~ x\), and 
f2 = 1 — fi- Otherwise, in those ranges of x for which the 
equilibrium state is single phase, it is possible, if desired, 
to perform a Legendre transform, so that the coefficients 
in the Landau free energy (a, /3, etc.) are expressed as 
functions of x and T. 




FIG. 1. Schematic phase diagrams derived from the Lan- 
dau free energy under the various conditions described in the 
text. "IAF" and "SC" indicate incommensurate (striped) an- 
tiferromagnetic and superconducting order, respectively. The 
circles represent classical critical or multicritical points and 
the squares quantum critical points. The various vertical lines 
represent trajectories through the phase diagram discussed in 
the text. The pale blue phase boundary in c) represents the 
effect of an applied magnetic field on the phase diagram. 

Generally, broken symmetry phases occur only at lower 
temperatures, so it is reasonable to expect a and a to 
change sign at bare transition temperatures defined ac- 
cording to a(x,T) = ao(x,T)[T — T s (x)] and a(x,T) = 
ao(x, T)[T — Ta(x)] where c*o, clq > 0. (Here, we use the 
same symbol for the original and Legendre transformed 
coefficients in the Landau theory.) On both theoretical 
and empirical grounds we expect the native magnetic or- 
dering temperature T s (x) to be a generally decreasing 



function of x, reflecting the frustration of hole-motion 
in the magnetically ordered state. However, a peak in 
T s (x) at a special "commensurate" value of x can occur 
when the period of the stripe order is a small integer 
times the underlying crystalline lattice constant; such a 
commensurability effect gives rise to the anomalous 
stability of the long period magnetism at x = 1/8, as 
shown schematically in Fig. 1. Ta(x) is known empir- 
ically to be a non-monotonic function of x, reaching a 
peak around an optimal value of x ~ 0.15 and dropping 
slowly as x is increased or decreased. 

There are several generic phase diagrams derivable 
from such a Landau free energy. Most simply, a tetra- 
critical point can occur at (T,fi) — (T*,/i*) which leads 
to a phase diagram of the type shown in Fig. la. Here, 
in addition to the phases with either magnetic or su- 
perconducting order, there is a homogeneous intermedi- 
ate phase with bulk coexistence of the two orders. We 
have included in this and all subsequent frames of Fig.l 
a second superconducting phase ( which is permitted, 
but not required in the simplest Landau theory) at hole 
concentrations x < 1/8. This is motivated by the ob- 
servation that many of the high temperature supercon- 
ductors exhibit multiple humps in the superconducting 
transition temperature plotted against hole concentra- 
tion, most notably the YBa2Cu307_,5 (especially when 
lightly Zn doped |l6|]) and La 2 _ a; Sr :i; Cu04 related com- 
pounds. As mentioned above, we have also indicated a 
peak in the magnetic ordering temperature at x = 1/8. 
The corresponding anomalous supression of the super- 
conducting transition temperature at this point shown 
in the figure is a consequence [pi of the peak in the mag- 
netic ordering under the assumption that 7 > 0, even if 
Ta{x) is monotonic in this range of x. 

The tetracritical point obtains as long as solutions 
exist to the simultaneous equations a(n* ,T*) = 
and a(/j,*,T*) = which at the same time satisfy 
the inequalities /3(/i*,T*) > 0, b(n*,T*) > 0, and 
f3(n*,T*)b(n*,T*) > 4 7 (^*,T*) 2 . If, on the other 
hand, the last of these inequalities is violated, i.e. 
if V ,T*)b(n* ,T*) < 2 7 (^*,T*), the tetracritical 
point is replaced by a bicritical point. Not only does that 
mean that there is no phase with bulk coexistence of the 
two orders, it also means that below the bicritical point, 
there is a region of two-phase coexistence, where a hole- 
poor, magnetically ordered and a hole-rich, supercon- 
ducting state coexist, as shown in Fig. lb. In reality, this 
coexistence must not be taken literally. Because of the 
long-range Coulomb interaction between holes, macro- 
scopic phase separation is thermodynamically forbidden. 
Where macroscopic phase coexistence would occur in a 
neutral system, a form of Coulomb frustrated phase sep- 
aration JnJ is expected, leading to a state which is in- 
homogeneous on an intermediate length scale. This also 
means that the two, coexisting phases are in microscopic 
proximity to each other, and hence that a modicum of su- 
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perconducting order will be induced in the magnetic re- 
gions, via the proximity effect, and conversely |17|. Such 
competition is a recurring theme in this problem; there 
are empirical and theoretical reasons to believe [ p~8|j40| 
that the long period antiferromagnetic stripe order, itself, 
is at least in part a consequence of Coulomb frustrated 
electronic phase separation on a smaller length scale. 

A third possibility shown in Fig. lc occurs if there is a 
tricritical point, where a(^*,T*) = /3(fi*,T*) = 0, while 
all the other coefficients (including certain higher order 
terms, not discussed explicitly) remain positive. This 
leads to a phase diagram of the sort shown in Fig. lc. 
Here, superconductivity manifests itself below a phase- 
boundary that terminates on the edge of the two-phase 
region in a critical end point. 

For completeness, we present a fourth possible phase 
diagram topology, shown in Fig. Id, for which a more 
thorough analysis of the free energy function is necessary. 
Here, instead of a multicritical point, we consider the oc- 
currence of a simple critical point, below which phase 
separation occurs into hole rich and hole poor phases, 
neither of which is ordered. In this case, both the anti- 
ferromagnetic and the superconducting phase boundaries 
terminate at critical end points. 

Relation to Experiment in L^CuO^a: The 
thermal evolution of a given material should be associ- 
ated with a trajectory in one of the generic phase dia- 
grams in Fig. 1. In particular, we propose associating 
La2Cu04+a with the vertical dashed line in Fig. lc or 
Id. It is a special trajectory, in the sense that it is tuned 
to pass close to the critical end point, but this requires 
fine tuning of only one parameter. Below T c , the system 
forms an inhomogeneous mixture of a high density super- 
conducting and a low density antiferromagnetic phase. 
At T c , the sample is a single-phase superconductor, with 
a superconducting volume fraction /g c (T) = X—fMag(T) 
that shrinks at the expense of the antiferromagnet as T 
is reduced through the two-phase region. Since there 
are no critical effects on the shape of the phase bound- 
ary associated with a critical end point, just below T c , 
fuag — A(T C ~ T) + . . ., where A is determined by the 
slope of the phase boundary. Since the magnetic order- 
ing of the hole poor phase would set in at a temperature 
well above T c , the ordered moment M(T) in the antifer- 
romagnetic fraction is immediately large, and essentially 
temperature independent! The growth of the magnetism 
is associated more with the growth of the hole-poor frac- 
tion rather than with the rise of the order parameter 
within the hole-poor regions. 

Not only is this scenario consistent with the simulta- 
neous onset of superconductivity and magnetism, it also 
reconciles the neutron scattering and /x-SR data, which 
we reproduce in Fig. 2. Specifically, in a two-phase mix- 
ture in which only one phase is magnetic, the intensity 
/ of the Bragg scattering measured by neutrons is re- 
lated to fMag and M, measured in /x-SR, according to 



the relation 



I(T) = M 2 (T)f Mag (T). 



(3) 
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FIG. 2. Interplay between magnetic and superconducting 
order in La2Cu04+^. The open and closed circles represent, 
respectively the antiferromagnetic fraction, fMag{T), and the 
product M 2 (T)fM a g(T) from the tt-SR data in Ref. @ . The 

ren squares represent the neutron intensity I(T) from Ref. 
, scaled by a factor of 2.7. The solid line is a theoretical 
prediction for FMag(T), using the lever rule, for the verti- 
cal dashed trajectory in the tricitical phase diagram in Fig. 
lc, assuming the two-phase region is bounded by the curve 
T2 — 4To(xi — x)(x — X2)/(x2 — x\) 2 with parameters (dis- 
cussed in the text) T = 43. 7K, x\ = 0.125, x 2 = 0.188, and 
with a mean hole density x = 0.15, representative Q of stage 
IV La 2 Cu0 4 +i. 

The fact that the absolute neutron intensity Q at low 
temperature is 0.022 /i^, which is 40% of that antici- 
pated from the muon experiments according to Eq. |3| 
suggests that either the sample interior penetrated by 
the neutrons is different from the surface region probed 
by the muons, or that the integration of the magnetic 
signal in the neutron experiment could be incomplete. 
In the figure, we have scaled the neutron data so that it 
matches the inferred intensity from /x-SR data at T = 0; 
the (scaled) neutron data still falls somewhat short of 
M 2 {T)fMag(T) immediately below T c , possibly because 
the neutron data are peak intensities, rather than in- 
tegrals over the three dimensional (in reciprocal space) 
structures containing the net spectral weight responsi- 
ble for the muon data. More specifically, if there is any 
broadening in the peaks as T c is approached, the peak 
intensity will be reduced relative to the integral. The red 
line in the figure corresponds to the lever rule prediction 
for an average hole density ^ x = 0.15 and assuming 
the simplest possible parabolic form for the bounding 
curve of the 2-phase region T%{x) — 4Tq{x% — x)(x — 
Xi)/(x2 — xi) 2 . So as to have no free parameters in de- 
termining the theoretical curve, we have taken x± = 1/8, 
reflecting the special stability of the striped phases at 
x w 1/8. The remaining parameters can be deduced 



3 



directly from f = /mos(O) = 0.4 and the condition 
that antiferromagnetism and superconductivity onset at 
the same temperature, T 2 (x) = T c = 42K, according to 
T = T c /[4/o(l - Jo)] and x 2 - x t = {x - Xl )/f . It 
turns out that the theoretical curves are very sensitive 
to the exact value of Jq. For instance, the quality of the 
fit can be improved if we take /o = 0.35 which is some- 
what smaller than the reported value, although possibly 
within experimental uncertainty. 

In the discussion above, we have ignored non-electronic 
physics. The most obvious possibility here is that mo- 
tion of excess oxygen results in an inhomogeneous dis- 
tribution of oxygen, which in turn would lead to inho- 
mogeneous hole density and electronic properties. The 
final outcome would be an inhomogeneous distribution 
of magnetism and superconductivity, even while charge 
neutrality would obviously be satisfied on a local scale. 
Nonetheless, our sense is that this explanation of the 
data is improbable. First, the onsets of superconductiv- 
ity and magnetism coincide |^| in La2- a: Sr a: Cu04 as well 
as La2CuC>4 + «5 in the same range of average hole density. 
Second, the findings of Lee et al Q provide an impor- 
tant clue concerning how charge neutrality is preserved 
at long length scales without substantial oxygen motion; 
they report magnetic order with a remarkably long cor- 
relation length (greater than 125A) in the basal plane, 
but with interplanar correlations extending only over 2-3 
planes. Charge neutrality can therefore be preserved over 
distances of order 10-15 Angstroms, even while the sys- 
tem breaks up into thin magnetic and superconducting 
layers - "pancakes" parallel to the basal planes. Finally, 
NMR studies [^0| suggest that the material is electron- 
ically single-phase at temperatures above T c , implying 
that the observed inhomogeneities are induced by the 
onset of order. 

La2- a; Sr 2 ,Cu04, Lai^-icNdo^Sr^CuC^, etc.: The 
phase diagrams in Fig.l provide a framework for under- 
standing many other properties of the lanthanum cuprate 
family. To begin with, a miscibility gap leading to co- 
existence of superconducting and non-superconducting 
phases readily accounts for the finding of optimal Meiss- 
ner fractions for La2_ a; Sr 2 ;Cu04 only near special hole 
densities |21f| . In addition, as for the superoxygenated 
La2Cu04+a discussed above, the ordered magnetic mo- 
ments deduced from neutron diffraction B are less than 
the frozen local moments deduced from muon spin re- 
laxation Q and much less than that seen for ordi- 
nary insulating two-dimensional antiferromagnets, im- 
plying also that the magnetic order resides in only a 
part of the sample. The appropriate phase diagram for 
La2_ a; Sr a ;Cu04 might then look like a disorder broad- 
ened (glassy) image of Fig. lb. On the other hand, 
as Nd is inserted, the magnetism (as detected in neu- 
tron diffraction) becomes stronger |j]-|3| and seems to 
appear throughout the sample volumes H, even while 
superconductivity survives. There is also remarkable 



evidence for a non-monotonic temperature dependence 
of the superfluid density [^2), which implies that on 
cooling, Lai.55Ndo.3Sro.i5Cu04 first undergoes a tran- 
sition to a uniform superconducting state, and then to 
a state with coexistence of magnetism and supercon- 
ductivity. Thus, the tetracritical diagram, Fig. la might 
be more appropriate for Lai.6- a; Ndo.4Sr 2 ;Cu04, with a 
microscopic coexistence of magnetic and superconduct- 
ing order; Lai.5sNdo.3Sro.i5Cu04 might then be repre- 
sented by the solid brown trajectory in that figure. Fi- 
nally, Lai-aBa^CuC^ exhibits two separated supercon- 
ducting "domes" , with an intervening magnetic regime 
(at x = 1/8) which is magnetic and not superconducting 
fl23H26| ■ The corresponding phase diagram could there- 
fore be that shown in Fig. la or lb. Of course, any real 
material exists not on a one-dimensional axis represent- 
ing the doping, but rather in a multidimensional space 
spanned by the parameters required to shift from dia- 
gram to diagram in Fig. 1. The outcome is then that, 
as demonstrated by inelastic neutron scattering from 
Lai.86Sro.i4Cu04 ^7j, optimally doped La2- a: Sr a ;Cu04, 
described by the dotted trajectory in Fig la, can show be- 
havior associated with the quantum critical point where 
magnetism disappears in the tetratcritical diagram. 

Additional Details: With the exception of 
La2Cu04+5, all the materials discussed have intrinsic 
disorder due to the random arrangement of the dopant 
atoms; any stripe ordering transition is thus expected |3~i| ] 
and observed |55| to be intrinsically glassy, with the or- 
dered phase being a "stripe glass." It should also be clear 
that the interplay between magnetism and supercon- 
ductivity, which accounts for so many key observations, 
omits other features |?6]-[59| expected or observed, of 
the actual phase diagrams. Within the magnetic regime 
there can be a variety of phases which we have not indi- 
cated: the long period magnetic order can be commensu- 
rate with the underlying lattice or incommensurate, and 
the stripes can point along the copper oxide bonds ( "hor- 
izontal") or (28) at 45° to them ("diagonal"). Finally, 
we have not shown the transitions involving charge or- 
der although certainly, at least for Lai.6~ a; Ndo.4Sr ;E Cu04, 
there is a separate transition Jl],[| at which unidirectional 
charge density wave ("charge stripe") order appears at 
a temperature above the magnetic and superconducting 
transitions. 

Prospects: A miscibility gap in the phase dia- 
gram not only resolves old puzzles, but also provides 
a framework for understanding current and future ex- 
periments. Particularly important are measurements of 
magnetic field-dependent effects, as they permit a con- 
tinuous variation of the parameters in the Landau free 
energy. For example, if La2Cu04+a in zero field hap- 
pens to lie on a trajectory which passes through a criti- 
cal end point, as we have supposed, then in a magnetic 
field, which will suppress the superconducting T c , mag- 
netic Bragg scattering, originating from a small magnetic 
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frac tion, will still appear at a temperature roughly equal 
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to t jic zero field T c . This is illustrated in Fig, lc, where 



mat/0007158. See also Weidinger, A., Niedermayer, Ch., 



the field shrinks the superconducting region of the phase 
diagram from the solid to the lighter blue line. How- 
ever, as we pointed out above, the vortex state is compli- 
cated - it is clear that physics beyond the simple Landau 
theory needs to be invoked |2^] in order to understand 
the dramatic magnetic field-induced increases in the an- 
tiferromagnetic Bragg intensities. This physics is clearly 
beyond the scope of the present work, but may relate 
to early theory indicating that vortices in superconduc- 
tors derived from Mott-Hubbard insulators are insulating 
nano-antiferromagnets, with a different charge density 
than the surrounding superconductor [j30| . Magnetore- 
sistance data |nj indicate enhanced insulating tenden- 
cies for x near 1/8, i.e. fields above H C 2 uncover the 
behavior also seen when superconductivity is suppressed 
by chemical pressure. Furthermore, neutron diffraction 
reveals field-induced magnetic Bragg scattering, which 
sets in near the zero-field critical temperature for super- 
conductivity f33f| . Even for samples beyond the miscibil- 
ity gap, because type II superconductors below H C 2 are 
heterogeneous mixtures of "normal" vortices and super- 
conducting material, the vortices can exhibit - on a finite 
length scale - the magnetism one might have expected if 
superconductivity had not intervened (assuming 7 > 0). 
Recent inelastic neutron scattering experiments jlij on 
optimally doped I^-ajSr^CuCU are in agreement with 
this expectation - the vortices are found to behave as 
nanomagnets with growing "stripe" order with decreas- 
ing temperature. 
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